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The discovery of the cosmic acceleration [J] has stim¬ 
ulated the search for models of the so called dark energy 
J2] responsible for this phenomenon. The crucial feature 
of the dark energy is that w = p/e < —1/3, where p 
is the pressure and e is the energy density. Some ob¬ 
servations p] hint to the possibility that the equation of 
state parameter w < — 1. The corresponding models are 
called phantom dark energy ones 0. These models have 
some unusual properties: to realize them one often uses 
the phantom scalar field with the negative sign of kinetic 
term; in many models the presence of the phantom dark 
energy implies the existence of the future Big Rip cos¬ 
mological singularity ||; according to some observations 
the crossing of the phantom divide line w = —1 occurs, 
the theoretical explanation of this fact also presents some 
kind of challenge J§]. 

The phantom model building has involved many differ¬ 
ent ideas. Here we would like to present a rather simple 
and natural cosmological toy model, linked to and in¬ 
spired by such an intensively developing branch of quan¬ 
tum mechanics and quantum Held theory as the study 
of non-Hermitian, but CPT (or PT) symmetric models 
(7, 8(,{|, lOj. The main point of this approach consists in 
the fact that there exists a large class of non-Hermitial 
Hamiltonians, which nevertheless possesses real and of¬ 
ten positive definite spectrum. As found by Bender and 
Boettcher in their seminal paper [J], they are character¬ 
ized by a potential which in one-dimensional case sat¬ 
isfies the property of PT - invariance V(x ) = V*{—x). 
Non-trivial generalization to quantum field theory has 
also been considered 0. It has been suggested that 
non-Hermitian quantum theory may find applications in 
quantum cosmology M- 

Here, we explore the use of a particular complex scalar 
field Lagrangian, which has real solutions of the classical 
equations of motion. Thereby we provide a cosmological 
model describing in a natural way an evolution from the 
Big Bang to the Big Rip involving the transition from 
normal matter to phantom matter, crossing smoothly the 


phantom divide line. The interest of our approach is 
related to its focusing on the intersection between two 
important fields of research, hopefully allowing for a their 
mutual cross-fertilization. 

In particular, we give an example of charged scalar 
matter interacting with a non-Hermitian potential which 
however does not break the CPT symmetry. In our 
model the classical solutions in the presence of gravity 
( FRW cosmological background) are such that the orig¬ 
inally complex Lagrangian becomes real on classical vac¬ 
uum configurations while one of the scalar component 
obtains the ghost sign of kinetic energy. Thereby we re¬ 
cover a more conventional phantom matter starting from 
the complex matter with normal kinetic energy. 

Let us consider a CPT symmetric, but non-Hermitian 
Lagrangian of a scalar field 

L = W - (1) 

with a potential V((f>, </>*) satisfying the CPT invariance 
condition 

(F(^0*))* = V(0*,0), (2) 

while the condition 

(n/>,<nr = m<n, (3) 

is not satisfied. (Indeed, it should be noted that P- 
invariance is trivial for a scalar field). For example, such 
potential can have a form 

= Fi(^ + <nw -p)- (4) 

If one defines 

= 1 + ifo (5) 

and considers potentials of the form 

V {$,(!>*) = bo(</>i) exp(*a</ 2 )) 


( 6 ) 
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where a is real parameter a, one can recognize the link 
to the so called PT symmetric potentials. 

Here, the functions <pi and (f> 2 are introduced as the 
real and the imaginary parts of the complex scalar field 
</>, however, in what follows, we shall treat them as inde¬ 
pendent spatially homogeneous variables depending only 
on the time parameter t. The equations of motion for 
fields 4>\ and (f> 2 have the form 

4>i + 3h<f>i + exp {iafo) = 0, (7) 


i(j> 2 + 3ih(f>2 - aV 0 ((j)i) exp(ia</> 2 ) = 0, (8) 


We provide now a simple realization of this idea by an 
exactly solvable cosmological model by implementing the 
technique for construction of potentials for a given cos¬ 
mological evolution It is convenient to start with a 
cosmological evolution as given by the following expres¬ 
sion for the Hubble variable: 


h{i) = 


A 

t{t R - t)' 


(17) 


The evolution begins at t = 0, which represents a stan¬ 
dard initial Big Bang cosmological singularity, and comes 
to an end in the Big Rip type singularity at t = t R . The 
derivative of the Hubble variable 


where h = - is the Hubble variable for a flat spatially 
homegeneous metric 

ds 2 = dt 2 — a 2 (t)dl 2 , (9) 

satisfying the Friedmann equation 

h 2 = \<t>\ + + v o((t>i) exp(za<^2)- (10) 

Let us notice, that the system of equations 171). 151. (fTUl) 
can have a solution where 4>i(t) is real, while the <f >2 is 
imaginary, or, in other words 


vanishes at 


• _ A(2t - t R ) 
h ~ t 2 (t R - t) 2 



(18) 


(19) 


when the universe crosses the phantom divide line. 

Next, we can write down the standard formulae con¬ 
necting the energy density and the pressure to the Hubble 
variable and its time derivative: 


Mt) = -m, (n) 

where £(f) is a real function. In terms of these two real 
functions, our system of equations can be rewritten as 


2 


£ 

2 


V o (0i)e°* = h 2 


A 2 

t 2 (t R ~t) 2 ’ 


( 20 ) 


^i+ 3 \/^i “ +K)(0i)exp(aC)^i+V o , ( < ? !, i) ex P( Q! C) = 0, 

( 12 ) 


f-V-Wi )S* = -±h-h 2 = - 


A(4t - 2 t R + 3 A) 
3 t 2 (t R - t) 2 


The expression for the potential Vq {<Pi ) follows 


( 21 ) 


- t }£ 2 + Vb(^i)exp(a^)^-aHo(0i)exp(a^) = 0. 

(13) 

Now, substituting </> 2 (t) from Eq. ITT1) into Eq. El 
we have the following expression for the energy density 

£ = h 2 = ^4>\ - ^i 2 + V 0 ((f>i)exp(a€). (14) 

The pressure will be equal 

p=\<t> 1 - \i 2 - v o{<t>i) exp«)- (15) 

It is easy to see that if < £ 2 the pressure will be 
negative and p/e < — 1, satisfying the phantom equation 
of state. Instead, when (j)\ > £ 2 , the ratio between the 
pressure and energy density exceeds —1 and, hence, the 
condition 

# = ^ (16) 

corresponds exactly to the phantom divide line, which 
can be crossed dynamically during the evolution of the 
field components and £(f). 


Wi) 


A(2t — t R + 3H) 

3 t 2 (t R -ty e 


The kinetic term satisfies the equation 


2 ^2 _ 2A(2t - t R ) 
11 * 3t 2 (t R ~t) 2 ' 


( 22 ) 


(23) 


It is convenient to begin the construction with the so¬ 
lution for £. Taking into account the formulae 0 and 
(1221) Eq. 1131 1 can be rewritten as 


I + 3C 


A 

t(t R - 1) 


aA(2t — t R T 3H) 
3t 2 {t R -1 ) 2 


(24) 


Introducing a new parameter 



(25) 


Eq. |2l looks like 


mt R 

v+v W^t) 


amt R (2t + t R (m 



( 26 ) 


Qt 2 (t R -t) 2 
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where 

y = i (27) 

It is not difficult to show that the solution of Eq. m is 
given by 


For the case a 2 m = 18 the function cpi(t) can be easily 
found from Eq. and it looks like follows: 


= ±\/32Arctanh 


tR - t 


tR 


(37) 


amt R (tR - t) m 

9 t m 


dt 


(2 1 + (m - 1 )t R )t m ~ 2 
(t R - t) m+2 


(28) 


where the inessential constant of integration will be dis¬ 
regarded. 

Let us estimate the behavior of the solution OBJ at 
t —> tR. Simple estimation gives 


am. 

9 (t-R - t)' 


(29) 


On the other hand the equality EH should be satisfied 
for all the values t < tR. That means that the value of 
£ 2 should be grater than the absolute value of right-hand 
side of Eq. ED- This last quantity at the limit t —* tR 
behaves as 2m/9(tR — t) 2 . Thus, one should have 


a 2 m 2 2m 

81 (t R - t) 2 ~ 9 (t R - t) 2 


(30) 


or, in other words, 


m > (31) 

or 

Before considering the concrete values of to, notice that 
the equation of state parameter w in the vicinity of the 
initial Big Bang singularity behaves as 

w = -l + —, (32) 

m 

while approaching the final Big Rip singularity this pa¬ 
rameter behaves as 


w = — 1 ——. (33) 

m 

Notice that the range for w does not depend on a, de¬ 
pending only on the value of the parameter to, which 
relates the scales of the Hubble variable h and of the 
time of existence of the universe tR. 

Remarkably, an integral in the right-hand side of Eq. 
EH is calculable analytically 


• amtR 
9 t(t R -1) 


(34) 


while 

fyrrt 

£=-g-(log*-log(ti{-t)). (35) 


From now on the parameter t will be dimensionless. In¬ 
clusion of characteristic time does not change the struc¬ 
ture of the potential because of its exponential depen¬ 
dence on £. Substituting the expression m into Eq. 
EH one has 


mtR((a 2 m. + 18 )tR — 36 1) 


(36) 


One can choose the positive sign in Eq. G3 without 
loosing the generality. 

Inverting Eq. m we obtain the dependence of the 
time parameter as a function of <f>i 


t = 


tR 


r osh 2 -4i= 
cos 


(38) 


Substituting expressions 03) and 03 into Eq. 
we can obtain the explicit expression for the potential 


Wi) 


2cosh»^j(2 +1 7 C c,sh 2 ^) 


t 


2 

R 


(39) 


We would like to emphasize that this potential is real 
and even. It is interesting that the time dependence of 
<f>i (t) could be found also for an arbitrary value of the 
parameter to, but for a 2 m > 18 this dependence cannot 
be reversed analytically and, hence, one cannot obtain 
the explicit form of the potential Vo(<j>i). 

Now, let us turn to Eq. 03, expressing the depen¬ 
dence of the time parameter t on <pi. It is convenient 
to consider the evolution of the value of <j >i between the 
values — oo, corresponding to the initial cosmological sin¬ 
gularity at t = 0, and 0, corresponding to the Big Rip 
singularity at t = t^. At the moment tp = tR/ 2, the 
equality £ 2 = (j)\ is satisfied and the universe is crossing 
the phantom divide line. It is easy to obtain from Eqs. 
03 and 03 the values of the fields £ and <j>\ at t = tp: 


£,(t P ) = 0, 


4>\{tp) = -\732Arctanh 


1 

2 ' 


(40) 


The potential m is smooth together with all its deriva¬ 
tives at this point. Moreover, it is quite regular at all the 
finite values of the scalar field 4>\. It is curious that the 
potential Vo(</>i) is finite also at the moment of the Big 
Rip. It is not, however, strange, because it enters in the 
expressions for the energy density 03 and the pressure 
03 being multiplied by the factor e a ^ which is singular 
at t = tR. 

In conclusion we would like to stress that we have con¬ 
structed a model relaxing the requirement of Hermiticity 
of the Hamiltonian of the theory which is equivalent to 
the reality of the classical Lagrangian. This relaxation, 
however, does not imply the breakdown of Lorentz and 
CPT invariance. For our classical solutions, expressed in 
terms of real fields, observable quantities like energy den¬ 
sity, pressure, Hubble variable turn out to be real. As a 
consequence our model describes in a rather natural way 
the transition from normal matter to phantom one. 


8lt 2 (t R - t) 2 
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